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INTRODUCTION 

Breast,  prostate,  Ixmg  and  renal  caneers  are  the  most  common  primary  tumors  which 
metastasize  to  bone.  The  vertebral  column  is  the  most  frequent  skeletal  site  for  metastases  (29)  and 
due  to  the  proximity  to  the  spinal  cord,  5%  to  10%  of  all  cancer  patients  develop  neurologic 
manifestations  (3).  Up  to  1/3  of  all  cancer  patients  develop  metastases  to  the  spinal  column  and 
over  50%  of  spinal  metastases  with  neurologic  manifestations  in  females  are  found  to  arise  from 
primary  breast  neoplasms  (3).  For  instance,  in  a  15  month  study  conducted  by  Toma  et  al.  (28) 
74.5%  of  patients  with  malignant  tumors  and  secondary  bone  involvement  had  a  primary  tumor 
located  in  the  breast,  of  which  54.5%  had  secondary  spinal  involvement.  Unfortunately,  there  is 
no  established  method  for  assessing  the  prognosis  of  patients  with  breast  cancer  and  metastasis 
confined  initially  to  bone  (30),  yet  bone  metastases  are  a  major  source  of  morbidity  in  patients 
with  breast  cancer  (i3).  In  the  face  of  greater  life  expectancy  for  women  with  breast  cancer,  spinal 
metastases  are  an  increasing  source  of  morbidity  which  negatively  affects  patients’  long  term 
quality  of  life. 

Yearly  there  are  about  18,000  new  cases  of  metastatic  involvement  of  the  vertebral  body 
in  the  United  States,  often  resulting  in  vertebral  instability  (i).  Metastases  to  the  spinal  column 
may  be  detected  by  radiographs,  bone  scan  or  may  present  with  symptoms  of  pain,  paralysis, 
loss  of  sensation  or  bowel  and  bladder  control.  Fracture  risk  prediction  has  significant  clinical 
importance,  as  prevention  of  fi’acture  in  high  risk  patients  may  be  possible  through  medical 
prophylaxis,  use  of  external  bracing,  or  internal  stabilization.  While  the  majority  of  vertebral 
body  fractures  occur  anteriorly,  posterior  body  (burst)  fractures  can  arise  from  high  impact 
loading  or  normal  loading  in  bones  compromised  by  osteoporosis  or  tumor.  Burst  fracture  is  one 
of  the  most  difficult  injuries  of  the  spine  to  successfiilly  treat,  in  part  because  the  exact 
mechanism  by  which  the  distraction  forces  are  transmitted  to  the  intracanal  fragments  of  the 
burst  fracture  have  not  been  adequately  investigated  (8).  The  average  survival  time  of  breast  cancer 
patients  with  vertebral  involvement  is  99.4  months,  or  40.4  months  from  the  time  of  diagnosis  of 
the  bone  metastases  (2).  Longer  survival  times  for  patients  with  metastatic  disease,  especially 
from  primary  breast  tumors,  have  increased  the  incidence  of  pathologic  fracture. 

While  improved  imaging  modalities  allow  clinicians  to  better  determine  the  extent  of  bony 
metastases,  assessing  the  risk  of  neurologic  compromise  is  guesswork  at  best.  In  a  review  of  70 
cases  of  spinal  cord  compression  in  breast  cancer.  Hill  et  al.  (9)  found  that  all  patients  had 
radiologic  evidence  of  bone  metastases  at  the  time  of  spinal  cord  compression.  The  majority  of 
their  patients  had  warning  symptoms  of  spinal  cord  compression  and  nearly  all  had  evidence  of 
vertebral  metastases  before  compression  occurred.  Accurate  neurologic  risk  assessment  would 
benefit  patients  by  providing  more  aggressive  treatment  options.  Clinical  and  radiographic 
studies  have  concluded  that  early  and  aggressive  treatment  of  vertebral  body  metastases  should 
be  considered  to  improve  patients’  life  quality  (9,28).  Elective  stabilization  with  posterior 
instrumentation  and  fusion  in  high  risk  patients  is  preferable  to  emergent  decompression  and 
uncertain  neurologic  recovery.  In  patients  with  burst  fractures  resulting  in  spinal  cord  or  cauda 
equina  compression,  complete  recovery  following  surgical  intervention  occurs  in  only  1/3  of 
patients  (26).  Although  neurologic  recovery  was  better  in  younger  female  patients  suffering  from 
metastatic  breast  cancer  than  in  older  men  with  prostatic  metastases  (26),  prevention  of  vertebral 
body  fracture  and  spinal  cord  compression  would  considerably  reduce  the  morbidity  arising  from 
vertebral  body  metastases.  Especially  in  cases  of  long  life  expectancy  with  tumors  like  breast 
cancer,  posterior  decompression  and  stabilization  can  exert  a  long-term  therapeutic  effect  (22).  A 
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better  understanding  of  spinal  canal  risk  may  permit  some  patients  to  avoid  surgery  altogether,  an 
important  goal  in  patients  with  limited  life  expectancy. 

Investigations  have  looked  at  the  relative  contributions  to  load  sharing  as  well  as  age 
related  changes  and  other  pathologic  processes  affecting  the  spine,  yet  most  have  focused  on 
anterior  failures  rather  than  posterior  body  (burst)  fracture  patterns.  Previous  methods  used  for 
strength  assessment  in  the  vertebral  body  have  included  finite  element  modeling,  computed 
tomography  and  in  vitro  mechanical  testing  (6,12,16,18,23.24,27).  The  non-invasive  techniques  can  be 
employed  to  provide  clinical  information  on  the  location  and  magnitudes  of  defects  and  stresses 
in  bones  and  for  prediction  of  fracture  risk  (10.17).  Determination  of  vertebral  body  strength  can 
provide  information  on  the  risk  of  vertebral  body  failure  including  fracture  location  and  potential 
for  canal  encroachment.  Previous  work  has  identified  metastatic  breast  cancer  lesions  at  risk  of 
fracture  in  the  femur  with  retrospective  clinical  record  and  radiographic  examination  (is)  and  for 
lesions  in  the  femoral  neck  with  finite  element  analysis  validated  through  mechanical  testing  (s). 

Previous  studies  have  attempted  to  predict  fracture  risk  of  metastatically  involved 
vertebral  bodies  using  defects  modeled  as  voids,  neglecting  any  effects  from  the  mass  and  fluid 
behaviour  of  tumor  material  (6,19,20,25).  Mizrahi  et  al.  (20)  developed  a  finite  element  model  to 
study  the  effects  of  metastatic  lesions  on  the  vertebral  body  strength.  The  study  used  a 
symmetric  model  and  represented  the  defects  simply  as  voids.  They  found,  not  surprisingly,  that 
larger  defects,  particularly  when  involving  greater  than  40%  of  the  vertebra,  were  likely  to  result 
in  greater  risk  of  fracture  in  the  anterior-superior  part  of  the  vertebral  body.  The  study,  however, 
did  not  address  the  effects  these  defects  have  on  creating  a  potential  risk  of  burst  fracture  and 
associated  canal  compromise. 

An  in  vitro  experimental  study  of  metastatic  lesions  in  the  thoracic  vertebrae  was 
performed  by  McGowan  et  al.  (i9)  to  predict  that  strength  reductions  were  associated  with  a 
reduction  in  vertebral  cross  sectional  area.  Mechanical  testing  of  vertebral  bodies  with 
experimentally  produced  trabecular  bone  defects  was  performed  under  axial-flexion  loads.  The 
results  suggested  that  strength  reduction  due  to  lytic  defects,  considered  as  voids  within  the 
thoracic  vertebrae,  is  proportional  to  the  cross  sectional  area  of  the  removed  bone.  However,  in 
no  case  was  the  posterior  vertebral  body  wall  found  to  be  displaced  into  the  canal  after  failure. 
Similar  mechanical  testing  was  performed  in  vitro  by  Silva  et  al.  (25).  They  found  significantly 
lower  failure  loads  for  vertebrae  with  lesions  regardless  of  the  position  of  the  defect.  However  all 
failures  occurred  in  the  anterior  vertebral  body  without  damage  of  the  posterior  wall,  thus  again 
the  true  effect  of  the  defects  on  posterior  stability  and  the  risk  of  neurologic  deficit  was  not 
addressed.  Dimar  et  al.  (6)  also  recently  tested  spinal  motion  segments  containing  defects  drilled 
as  a  cylinder  through  the  trabecular  centrum  as  well  as  the  anterior  and  posterior  cortical  walls. 
This  removal  of  posterior  cortex  eliminates  the  possibility  of  assessing  the  risk  of  posterior 
fracture  or  neurologic  compromise. 

In  Mizrahi’s  finite  element  model  and  the  mechanical  testing  of  McGowan  et  al.,  Silva  et 
al.  and  Dimar  et  al.,  the  defects  were  modeled  as  voids,  neglecting  any  effects  which  may  arise 
due  to  the  mass  and  fluid  behaviour  of  tumor  material.  Clinically,  encroachment  of  either  tumor 
material  or  displaced  bone  may  result  in  canal  compromise,  thus  the  mass  effect  of  the  tumor  is 
likely  to  be  significant.  Consequently,  although  these  studies  give  some  insight  into  the  type  of 
lesion  which  may  put  a  patient  at  risk  for  compression  fracture,  they  do  not  consider  which 
lesions  pose  a  risk  for  neurologic  deficit. 
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The  mechanism  of  burst  fracture  has  been  clinically  suggested  to  be  due  to  entry  of  the 
nucleus  of  the  disc  into  the  vertebral  body  through  an  endplate  fracture  causing  the  body  to 
pressurize  and  ultimately  burst  (ii),  thus  indicating  poroelastic  behaviour  (the  rate  dependent 
behaviour  of  a  tissue  containing  both  a  solid  and  a  fluid  phase)  .The  two  phase  porous  behaviour 
of  trabecular  bone  has  been  demonstrated  by  Carter  and  Hayes  (4)  through  confined  compression 
testing  of  bone  with  and  without  marrow.  The  presence  of  marrow  increased  the  strength, 
modulus  and  energy  absorption  of  the  specimens  when  tested  at  a  strain  rate  which  would 
simulate  loading  conditions  found  in  high  energy  burst  fractures.  Downey  et  al.  (?)  found  an 
increase  in  intraosseus  pressure  in  an  intact  bone  structure  under  compressive  loading.  Shirado  et 
al.  (23)  tested  vertebral  body  motion  segments  (vertebra-disc-vertebra)  under  axial  compression 
and  were  able  to  reproduce  such  burst  fracture  patterns.  These  studies  examined  solid  specimens 
of  trabecular  bone  and  vertebral  bodies  devoid  of  any  defects.  In  metastatically  involved 
vertebral  bodies  however  the  tumor  mass  consists  of  a  solid  matrix  with  greatly  reduced  elastic 
stiffness  and  a  large  percentage  of  fluid.  This  large  fluid  mass  may  cause  the  effect  of  loading  rate 
on  bone  fracture  to  be  more  pronounced.  In  vertebral  bodies  compromised  by  tumor,  burst 
fracture  patterns  are  found  to  occur  under  physiologic  loading  conditions  similar  to  that  which 
occur  in  healthy  vertebrae  only  under  high  energy  impact.  Yet  an  analysis  of  the  poroelasticity 
of  trabecular  bone  within  intact  vertebral  bodies,  with  or  without  metastatic  involvement,  has  not 
been  studied.  Nor  have  any  studies  on  metastatically  involved  vertebral  bodies  included  the  two 
phase  nature  of  the  tumor.  If  metastatically  involved  vertebral  bodies  exhibit  poroelastic 
behaviour  under  physiologic  loading  rates,  fracture  prediction  estimates  from  earlier  models  may 
not  accurately  represent  the  stress  distributions  in  the  vertebral  body. 

Guidelines  are  needed  to  permit  clinicians  to  select  appropriate  patients  for  intervention 
before  the  development  of  neurologic  compromise.  Current  guidelines  for  prophylactic 
stabilization  have  been  extrapolated  from  clinical  experience  and  subjective  intuition  rather  than 
from  experimental  evidence.  Suggested  indications  for  stabilization  include:  >50%  vertebral  body 
destruction,  pedicle  destruction,  involvement  of  two  or  more  adjacent  vertebral  bodies, 
involvement  of  anterior  and  posterior  elements,  progressive  collapse  of  a  vertebral  body,  and 
severe  pain  not  relieved  by  other  treatments.  A  widely  accepted  definitive  set  of  guidelines  for 
prophylactic  intervention  in  metastatically  involved  vertebral  bodies  has  yet  to  be  determined. 

The  purpose  of  our  research  is  to  quantify  fracture  risk  in  vertebral  bodies  affected  by 
metastatic  breast  cancer  in  order  to  develop  a  set  of  clinical  guidelines  for  the  prophylactic 
treatment  of  patients  with  such  metastases.  We  hypothesize  that  the  use  of  poroelastic  theory 
in  modeling  will  yield  more  realistic  estimates  of  fracture  prediction  compared  to  values  based  on 
a  linear  elastic  model  and  that  poroelastic  properties  of  both  bone  and  tumor  tissue  are  necessary 
to  understand  the  mechanisms  of  vertebral  body  failure  and  adequately  assess  the  risk  of 
neurologic  compromise.  The  ultimate  goal  of  this  research  is  to  use  engineering  analysis  to 
develop  clinical  diagnostic  and  surgical  guidelines  for  prophylactic  stabilization  of  metastatically 
involved  vertebral  bodies  in  order  to  prevent  paralysis  and  significantly  improve  the  quality  of 
life  of  patients  with  vertebral  body  metastases. 

In  our  initial  proposal  during  the  first  year  we  planned  to  determine  the  poroelastic 
properties  of  breast  tumor  tissue  and  develop  both  a  two  dimensional  axisymmetric  and  a  three 
dimensional  metastatically  involved  vertebral  body  model.  The  following  is  a  summary  of  our 
progress  towards  those  goals. 
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Technical  Objective  1:  Poroelastic  properties  of  tumor  tissue 

To  determine  the  poroelastic  properties  of  breast  tumor  tissue.  Tissue  properties  will  be 
measured  experimentally  through  confined  compression  testing  of  human  tumor  specimens 
obtained  at  surgery  from  bone  metastases. 

To  establish  the  mechanical  properties  of  the  metastatic  tumor  material  in  terms  of  linear 
poroelastic  theory  we  have  set  up  creep  experiments  on  tumor  specimens  in  confined 
compression  (Figure  1).  The  outcome  variables  for  each  specimen  are  the  aggregate  compressive 
modulus  (Ha)  and  the  hydraulic  permeability  (k).  Tumor  tissue  will  be  harvested  at  surgery 
from  patients  with  metastatic  involvement  of  bone,  and  frozen.  Cylindrical  test  specimens  (1-2 
mm  thick  x  5  mm  diameter)  are  prepared  by  first  cutting  out  blocks  of  1-2  mm  thickness  using  a 
vibratome,  and  then  punching  out  the  5  mm  diameter  plugs  directly  into  the  confining  chamber 
(Figure  2). 

Specimens  inside  the  cylindrical  confining  chamber  are  sandwiched  between  rigid,  porous 
glass  platens.  Normal  saline  is  added  to  the  chamber  and  the  compressive  creep  tests  performed 
in  the  bath  held  at  body  temperature  (37°  C).  A  constant  load  of  0.4  to  1  g  is  applied  to  the 
specimen  to  achieve  a  final  strain  of  <20%.  Tests  are  conducted  imtil  equilibrium  is  achieved  (1 
to  24  hours).  Due  to  possible  inaccuracies  in  the  load  application  through  the  counterweighted 
arm,  an  additional  load  cell  is  included  to  record  the  actual  force  felt  by  the  specimen.  The  creep 
displacement  versus  time  data  is  filtered  and  collected  using  LAB  VIEW  5.0  and  fit  to  a  solution 
derived  using  linear  poroelastic  theory  (necessitating  the  small  strain  requirements  of  the 
experimental  protocol).  The  parameters  Ha  and  k  are  determined  using  non-linear  best-fit 
regression  (MATLAB,  v5.0).  Based  on  clinical  availability  thus  far,  eight  tumor  samples  have 
been  collected  from  a  variety  of  tumor  metastases  to  bone.  Although  acquiring  a  sufficient 
number  of  specimens  for  the  study  has  taken  longer  than  expected,  we  do  not  anticipate  this  to 
be  a  limitation  in  the  final  project.  The  histology  of  each  tumor  is  analyzed  to  determine  the 
cellular  content  and  percentage  of  necrosis.  Results  of  stiffness  and  permeability  will  ultimately 
be  correlated  with  these  parameters  to  discover  whether  such  parameters  (available  from 
biopsy)  can  be  reliable  indicators  of  tumor  tissue  material  properties. 

Difficulties  in  the  execution  of  the  testing  protocol  have  occurred  due  to  the  very  small 
loads  needed  to  maintain  the  small  strain  conditions  necessary  to  implement  linear  poroelastic 
theory.  Initially  a  confined  compression  testing  set  up  was  designed  using  an  MTS.  However,  it 
became  evident  that  the  feedback  control  on  such  a  large  machine  was  not  adequate  to  maintain 
constant  loads  at  these  small  magnitudes.  As  such,  our  experimental  apparatus  was  designed  and 
built  to  apply  constant  load  without  the  necessity  of  a  variable  feedback  loop.  Further 
difficulties  have  resulted  due  to  friction  between  the  specimen  or  platens  and  the  confining 
chamber.  Due  to  the  small  applied  loads,  overcoming  such  frictional  forces  can  prove 
troublesome 

The  results  and  subsequent  curve  fit  to  determine  the  material  parameters  Ha  and  k  from 
one  specimen  of  metastatic  breast  tumor  tissue  are  shown  in  Figure  3.  While  only  a  very  limited 
number  of  specimens  have  been  tested,  thus  far  the  aggregate  moduli  of  the  tumor  tissue  tested 
have  been  within  the  range  of  0.001  to  0.005  MPa.  This  includes  both  specimens  from  primary 
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lung  and  breast  neoplasms.  However,  the  hydraulic  permeability  of  the  specimens  thus  far 
differ  by  several  orders  of  magnitude  from  values  in  the  range  of  10*’  to  10'^.  This  may 
ultimately  be  explained  by  the  difference  in  histology  of  the  specimens.  The  results  from  the 
breast  metastasis  tested  in  Figure  3  had  a  hydraulic  permeability  on  the  order  of  10"^  and 
consisted  of  a  highly  cellular  tissue  (up  to  75%  in  some  areas)  with  large  areas  of  necrosis.  In 
contrast,  one  lung  specimen  tested  demonstrated  a  hydraulic  permeability  on  the  order  of  10'^ 
and  was  made  up  of  a  considerable  amount  of  stroma  with  a  cellularity  of  approximately  30% 
and  less  than  10%  necrosis.  More  analyses  and  specimens  are  needed  before  and  conclusions  can 
be  drawn  about  the  material  properties  of  tumor  tissue  which  has  metastasized  to  bone. 


Technical  Objective  2:  Poroelastic  theory  applied  to  metastatically  involved  vertebral 
body  modeling 

To  develop  a  poroelastic  two-dimensional  axis-symmetric  finite  element  model  of  the  vertebral 
body  and  parametrically  assess  the  effects  of  rate  dependence  on  vertebral  body  strength  and  the 
implications  of  poroelastic  theory  in  the  consideration  of  metastatic  involvement  in  the  model. 

We  developed  a  two-dimensional  poroelastic  axisymmetric  finite  element  model  of  a 
spinal  motion  segment  consisting  of  the  first  lumbar  vertebral  body  and  adjacent  intervertebral 
disc  (Figure  4).  The  model  assumed  a  midsagittal  axis  of  symmetry  and  a  transverse  plane  of 
symmetry.  No  posterior  elements  were  modeled.  A  finite  element  mesh  consisting  of  493 
elements  (8-noded  quads  and  rebars)  was  generated  and  analyzed  using  commercial  software 
(PATRAN  5.0,  ABAQUS  5.6).  The  model  was  constructed  to  allow  the  inclusion  of  a  centrally 
located  tumor  in  the  vertebral  body,  representing  a  25%  or  50%  defect  of  the  trabecular  bone  by 
volume.  An  elastic  model  with  the  number  of  elements  increased  by  a  factor  of  four  was 
generated  and  demonstrated  mesh  convergence  was  reached  in  the  initial  model. 

The  intact  model  (no  tumor  included)  was  analyzed  under  a  fully  elastic  configuration 
(both  disc  and  vertebral  body  represented  with  elastic  material  properties),  with  a  poroelastic 
disc  attached  to  a  linear  elastic  vertebral  body  (with  and  without  restricted  flow  through  the 
cartilaginous  endplate),  and  as  a  fully  poroelastic  model.  In  the  analyses  of  the  metastatically 
involved  cases  only  poroelastic  modeling  was  performed. 

The  vertebral  body  was  composed  of  a  cortical  shell  and  endplate  modeled  isotropically 
and  a  transversely  isotropic  trabecular  bone  centrum,  with  properties  (both  elastic  and 
poroelastic)  chosen  to  represent  a  healthy  young  spine  (i4).  The  intervertebral  disc  consisted  of 
cartilaginous  endplate,  nucleus  pulposus  and  annulus  fibrosus.  The  nucleus  and  cartilaginous 
endplate  were  modeled  isotropically  in  the  elastic  and  poroelastic  configurations.  The  annulus 
was  composed  of  an  isotropic  material  modeled  hyperelastically  in  both  configurations  and 
reinforced  with  rebar  elements  representing  collagen  fibers.  The  tumor  tissue  was  modeled  as  a 
poroelastic  isotropic  material  with  the  aggregate  modulus  and  hydraulic  permeability  varied  to 
represent  a  spectrum  of  material  properties  formd  in  lytic  lesions.  In  the  poroelastic  and  mixed 
analyses,  the  fluid  phase  for  all  tissues  was  modeled  as  incompressible  with  properties  equivalent 
to  water. 

The  model  was  loaded  axially  under  a  uniform  pressure  of  1  MPa  applied  though  the 
midplane  of  the  intervertebral  disc  (which  corresponds  to  a  1200N  force,  approximating  the 
compressive  force  on  the  lumbar  spine  for  an  individual  standing  upright  holding  an  8.3  kg  mass 
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with  outstretched  arms  (21)).  The  intact  poroelastic  model  was  analyzed  under  a  range  of  loading 
rates  from  10,000  MPa/s  (impact)  to  4  MPa/s  (slow  walking).  In  the  mixed  elastic/poroelastic 
runs  a  physiologic  loading  rate  level  of  lOOMPa/s  was  applied.  In  the  analyses  of  the 
metastatically  involved  cases  poroelastic  modeling  was  performed  under  loading  rates  of  10 
MPa/s  and  lOOMPa/s. 

We  focused  on  trabecular  bone  failure  and  our  outcome  variables  were:  radial  displacement 
of  the  vertebral  body  midline  (Ul)  as  an  indicator  for  spinal  canal  encroachment,  endplate 
deformation  (U2)  as  a  measure  of  endplate  fracture,  maximum  strain  (E)  indicating  potential 
trabecular  bone  failure,  and  pore  pressure  (POR)  to  determine  the  amount  of  load  carried  by  the 
fluid  phase  in  the  poroelastic  analyses. 

When  we  represented  the  vertebral  body  elastically,  our  results  were  dependent  on  the 
disc  representation  used  (elastic  versus  poroelastic)  (Table  1).  Greater  disc  deformation  in  the 
axial  and  radial  directions  in  the  fully  elastic  model  create  different  loading  conditions  applied  to 
the  cortical  endplate  of  the  vertebral  body,  and  thus  yield  significant  differences  in  the  trabecular 
bone  outcome  parameters  (Figure  5).  However,  changes  in  the  boundary  conditions  at  the  disc  - 
vertebral  body  interface,  from  restricted  flow  to  free  flow  through  the  endplate,  did  not  result  in 
significant  differences  in  vertebral  body  response  in  the  mixed  models. 

The  fully  poroelastic  representation  yielded  higher  maximum  principal  tensile  strains  and 
radial  displacement  and  lower  maximum  principal  compressive  strains  and  axial  displacement  in 
the  vertebral  body  trabecular  bone  than  found  in  the  mixed  models.  Conversely,  results  within  the 
intervertebral  disc  under  the  fully  poroelastic  configuration  are  bounded  by  the  results  derived 
from  the  two  mixed  models. 

Tumor  size  had  the  greatest  effect  on  vertebral  body  displacement,  strain  and  pore 
pressure  in  the  poroelastic  model  (Table  2).  Increasing  tumor  size  from  25%  to  50%  yielded  a 
40%  increase  in  maximum  radial  displacement,  a  40%  increase  in  tensile  hoop  and  radial  strains 
located  at  the  axis  of  symmetry  along  the  bone  tumor  interface,  and  a  130%  increase  in  the 
compressive  strain  under  the  center  of  the  endplate.  Maximum  trabecular  bone  displacement  in 
the  radial  direction  shifted  from  an  area  adjacent  to  the  outer  endplate  in  the  intact  vertebral  body 
to  the  transverse  midline  when  a  defect  was  incorporated.  Pore  pressure  in  the  vertebral  body 
trabecular  bone  was  increased  by  160%  and  280%  respectively  with  the  inclusion  of  25%  and 
50%  tumors  into  the  vertebral  body  model  as  compared  to  the  case  with  no  tumor  involvement 
(Figure  6). 

In  the  intact  poroelastic  model,  increase  in  loading  rate  resulted  in  a  decrease  in  the 
maximum  displacement  in  the  radial  direction  at  the  distal  edge  of  the  model  adjacent  to  the 
endplate,  a  site  indicating  a  potential  increased  risk  of  compression  fracture  (Table  3).  This  was 
accompanied  by  an  increase  in  the  axial  displacement  under  the  endplate  adjacent  to  the  nucleus 
pulposus  for  loading  rates  up  to  0.01  s’^  This  axial  displacement  was  found  to  saturate  under 
loading  rates  faster  than  0.01  s'*.  As  well,  pore  pressure  within  the  vertebral  body  was  found  to 
decrease  with  increased  loading  rate  in  the  intact  model.  Conversely,  with  a  10  fold  reduction  in 
loading  rate  in  the  pathologic  models,  maximum  pore  pressures  in  the  vertebral  body  were 
reduced  by  8%  and  demonstrated  a  smaller  overall  pore  pressure  gradient  (Figure  4).  This 
reduction  in  loading  rate  in  the  pathologic  models  also  caused  maximum  endplate  displacement 
and  axial  compressive  strains  to  increase  by  5%  to  13%. 
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Variation  of  tumor  tissue  material  properties  did  not  have  much  influence  on  the  model 
outcome  parameters.  A  10  fold  increase  in  tumor  tissue  permeability  was  found  to  have  no  effect 
on  our  vertebral  body  results.  However,  a  25  fold  increase  in  tumor  aggregate  modulus  in  the  50% 
defect  model  did  result  in  small  reductions  in  vertebral  body  strain  (2%  to  6%),  displacement 
(0.5%  to  3%)  and  pore  pressure  (2%). 

In  assessing  the  need  for  utilizing  poroelastic  modeling  techniques  in  analyzing  spinal 
motion  segments,  it  was  found  that  the  boundary  conditions  applied  to  the  vertebral  body 
endplate  are  dependent  on  the  constitutive  assumptions  of  the  intervertebral  disc,  and  as  such, 
affect  the  strain  and  displacement  results  in  the  trabecular  bone  centrum.  Greater  disc 
deformation  in  the  axial  and  radial  directions  in  the  fully  elastic  model  create  different  loading 
conditions  applied  to  the  cortical  endplate  of  the  vertebral  body,  and  thus  different  responses 
within  the  vertebral  body  itself  For  analyses  aimed  at  studying  the  responses  of  the 
intervertebral  disc  alone,  our  results  suggest  it  is  reasonable  to  utilize  a  mixed  model.  That  is,  the 
results  found  in  the  intervertebral  disc  in  the  fully  poroelastic  model  are  bounded  by  the  results 
derived  from  the  mixed  models;  a  free  flow  boundary  at  the  endplate  simulating  a  fully 
poroelastic  model  at  a  slow  loading  rate,  and  a  no  flow  boundary  condition  simulating  fast 
loading.  However,  results  for  strains  and  displacements  in  the  poroelastic  vertebral  body  are  not 
bounded  by  the  mixed  model  results.  The  inclusion  of  the  fluid  phase  into  the  vertebral  body 
results  in  a  portion  of  the  spinal  load  being  supported  by  the  liquid,  reducing  the  solid  strains  in 
the  axial  direction  in  the  trabecular  bone.  Concurrently,  the  incompressibility  of  the  fluid  phase 
causes  expansion  and  higher  strains  to  develop  in  the  radial  direction  within  the  vertebral  body. 

In  determining  the  effect  of  loading  rate  on  a  spinal  motion  segment  it  was  found  that  as 
the  loading  rate  was  increased,  higher  values  of  compressive  strains  and  lower  values  of  tensile 
strains  were  predicted  in  the  solid  phase  of  the  vertebral  body  trabecular  bone.  At  the  same  time, 
maximum  pore  pressures  within  the  nucleus  were  found  to  increase,  while  pore  pressures  in  the 
vertebral  body  decreased.  This  opposing  trend  was  due  to  the  large  disparity  in  the  diffusion  time 
constants  of  the  nucleus  (on  the  order  of  10^  s)  versus  the  trabecular  bone  centrum  (on  the  order 
of  10'*  s).  As  fluid  does  not  have  time  to  escape  from  the  nucleus  imder  higher  loading  rates,  the 
incompressibility  of  the  fluid  phase  in  the  nucleus  causes  increased  axial  deformation  of  the 
vertebral  body  endplate.  These  results  agree  with  clinical  observations  indicating  an  increased 
likelihood  for  endplate  failure  in  vertebral  bodies  loaded  at  higher  rates,  potentially  leading  to 
burst  fracture  patterns  once  the  endplate  fracture  permits  the  nucleus  to  enter  and  pressurize  the 
vertebral  body. 

In  examining  the  relative  effects  of  tumor  size,  material  properties  and  loading  rate  in  the 
assessment  of  metastatically  involved  vertebral  bodies,  increased  tumor  size  was  foimd  to  cause 
the  greatest  increase  in  vertebral  body  displacements,  strains  and  pore  pressure.  The  location  of 
maximum  displacement  and  pore  pressures  suggest  a  greater  risk  for  endplate  and  radial  vertebral 
body  failure  (i.e.  an  increased  likelihood  of  vertebral  burst  fracture)  as  the  size  of  defect  increases. 
Greater  maximum  radial  displacement  values  and  increased  tensile  hoop  strains  predicted  at  the 
transverse  midline  of  the  vertebral  body  may  potentially  correspond  to  a  catastrophic  fracture 
pattern  with  bone  and  tumor  encroachment  into  the  spinal  canal  in  a  three-dimensional  model. 
This  finding  agrees  with  the  radiologic  assessments  of  Taneichi  et  al.  (Taneichi  et  ah,  1997),  who 
found  tumor  size  (along  with  pedicle  destruction)  to  be  the  most  important  risk  factor  for 
vertebral  body  collapse  in  metastatically  involved  spines.  Loading  rate  was  found  to  have  the 
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second  strongest  effect  on  metastatically  involved  vertebral  body  behavior.  With  inclusion  of  the 
poroelastic  defect,  increased  loading  rate  resulted  in  more  load  carried  by  the  fluid  phase  within 
the  centrum,  accounting  for  the  increased  pore  pressures  and  lower  axial  compressive  strains  and 
displacements  in  the  trabecular  bone.  Radial  displacement  along  the  transverse  midline  increased 
at  higher  loading  rates,  indicating  escalating  risk  of  potential  neurologic  compromise.  These 
findings  may  indicate  that  the  mechanism  of  burst  fracture  may  differ  in  intact  versus 
metastatically  involved  vertebral  bodies.  The  presence  of  a  defect  which  produces  increased 
vertebral  body  pressurization  and  radial  displacement  may  potentially  allow  a  burst  fracture 
pattern  to  occur  without  an  endplate  fracture.  Variation  of  tumor  tissue  material  properties  did 
not  have  a  large  effect  on  the  behavior  of  the  trabecular  bone  centrum  of  the  vertebral  body,  thus 
the  primary  site  of  the  lytic  tumor  may  not  be  of  major  significance  in  assessing  failure  risk. 

While  this  study  is  limited  by  the  use  of  a  two-dimensional  axisymmetric  model,  results 
of  this  simplified  analysis  provide  a  justification  for  the  use  of  more  complex  three-dimensional 
poroelastic  models  and  demonstrate  the  importance  of  loading  rate.  Due  to  the  use  of  an  idealized 
model,  the  results  found  here  are  used  solely  to  make  comparisons  between  the  different  cases 
(elastic/poroelastic,  loading  rate,  tumor  size,  tumor  tissue  material  properties).  Although  the  cost 
of  poroelastic  modeling  is  relatively  large,  both  in  terms  of  determination  of  additional  poroelastic 
material  properties  and  extended  computing  time  (especially  for  three-dimensional  models),  a 
purely  elastic  approach  may  not  be  sufficient  to  reproduce  the  complex  behavior  exhibited  by 
this  two  phase  structure. 

Task  3:  To  incorporate  poroelasticity  and  metastatic  defects  into  a  three-dimensional  vertebral 
body  model. 

We  have  developed  a  poroelastic  three-dimensional  model  of  a  spinal  motion  segment 
consisting  of  a  vertebral  body  and  adjoining  intervertebral  discs  cut  through  the  midsagital  plane 
(Figure  7).  The  model  includes  the  posterior  arch  but  no  other  posterior  elements.  It  has  been 
designed  to  incorporate  metastatic  defects  of  varied  size  (encompassing  15%,  30%  or  45%  of  the 
total  trabecular  volume,  and  with  or  without  involvement  of  the  pedicle)  in  order  to  allow  the 
determination  of  the  sensitivity  to  stress  and  posterior  vertebral  body  fracture  risk  of  varied 
configurations.  No  results  have  been  obtained  as  of  yet  from  analyses  of  the  model. 
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CONCLUSIONS 


Thus  far,  we  are  on  track  and  meeting  the  schedule  for  the  first  12  months  provided  in  our 
initial  proposal.  While  limited  numbers  of  specimens  have  been  tested  at  the  present  time,  the 
experimental  confined  compression  protocol  developed  appears  to  repeatably  measure  the 
poroelastic  material  parameters  of  tumor  tissue  specimens.  These  data  are  expected  to  correlate 
well  with  histologic  findings  in  comparing  tissue  cellularity  vs.  hydraulic  permeability  and 
aggregate  modulus.  This  information  should  prove  useful  not  only  to  contribute  to  our  further 
work  in  assessing  burst  fracture  risk,  but  also  in  other  areas  of  drug  delivery  and  cancer  treatment. 

Our  two-dimensional  finite  element  model  has  demonstrated  the  need  for  utilizing 
poroelastic  modeling  techniques  in  analyzing  spinal  motion  segments.  Faster  loading  rates  were 
shown  to  indicate  an  increased  likelihood  for  endplate  failure  in  intact  vertebral  bodies  in 
concurrence  with  clinical  observations.  Increased  tumor  size  was  found  to  cause  the  greatest 
increase  in  metastatically  involved  vertebral  body  displacements,  strains  and  pore  pressure  and 
suggested  a  greater  risk  for  endplate  and  radial  vertebral  body  failure  as  the  size  of  defect 
increased.  Higher  loading  rates  applied  to  the  metastatically  involved  models  also  indicated  an 
escalating  risk  of  a  potential  burst  fracture  pattern,  whereas  variation  of  tumor  tissue  material 
properties  did  not  have  a  large  effect  in  assessing  failure  risk.  These  simple  models  have  formed  a 
basis  which  has  guided  our  work  thus  far  in  three-dimensional  modeling  and  experimental 
validation  with  the  goal  of  determining  better  criteria  to  assess  burst  fracture  risk  in  patients  with 
spinal  metastases. 
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Table  1 :  Maximum  Displacement  (U)  and  Strain  (E)  Results 
in  the  Vertebral  Body  from  the  Elastic,  Mixed  Model  and 
Poroelastic  Analyses  (@  100  MPa/s). 


yi{mm) 

U2 

Ell  {x10'3) 

E22(xi0'3) 

E33  (xio*3) 

Elastic 

0.0493 

-0.0477 

3.56 

-4.91 

2.60 

Mixed  Free  Flow 

0.0479 

-0.0733 

3.58 

-6.14 

2.54 

Mixed  No  Flow 

0.0498 

-0.0720 

3.80 

-6.06 

2.64 

Poroelastic 

0.0561 

-0.0674 

3.78 

-5.63 

2.98 

Table  2:  Maximum  Displacement  (U)  in  the  Vertebral  Body  and 
Pore  Pressure  (POR)  Results  in  the  Vertebral  Body  and  Disc 
from  the  Poroelastic  Model  under  Varied  Loading  Rates. 


Loadino  Rate  fMPa/s^ 

yi{mm) 

U2  (mm) 

PORvb  (MPa) 

POROISC  (MPa; 

10000 

0.0488 

-0.0672 

0.294 

1.80 

1000 

0.0507 

-0.0673 

0.302 

1.80 

100 

0.0561 

-0.0674 

0.328 

1.80 

10 

0.0670 

-0.0627 

0.497 

1.80 

4 

0.0727 

-0.0580 

0.630 

1.77 

Table  3:  Displacement  (U,  mm),  strain  (E),  and  pore  pressure  (POR, 
MPa)  results  in  the  vertebral  trabecular  bone  of  intact,  25%  and  50% 
tumor  volume  models  at  loading  rates  of  10  and  100  MPa/s. 


Loadina  Rate 

Tumor  Volume 

Ul 

U2 

E11* 

E22* 

E33* 

POR 

100  MPa/s 

0% 

0.056 

-0.067 

0.378 

-0.563 

0.298 

0.328 

25% 

0.142 

-0.423 

2.21 

-1.35 

2.21 

0.832 

50% 

0.199 

-0.975 

3.08 

-0.878 

3.08 

1.24 

10  MPa/s 

0% 

0.067 

-0.063 

0.443 

-0.521 

0.358 

0.497 

25  % 

0.139 

-0.444 

2.23 

-1.43 

2.23 

0.768 

50% 

0.193 

-1.10 

3.50 

-0.996 

3.50 

1.15 

*  X10  -2 
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Figure  1:  Confined  compression  experimental  set-up 
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Figure  2:  Confined  compression  chamber  with  glass 
porous  platens,  LVDT  and  plunger. 
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Figure  4;  Finite  element  mesh  of  the  spinal  motion  segment 


Figure  5:  Deformed  mesh  of  the  fully  elastic  (left)  and 
poroelastic  (right)  representations  of  the  spinal  motion 
segment  under  a  1  MPa  load. 
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Figure  6  Pore  pressure  (POR)  results  developed  in  the 
vertebral  trabecular  bone  of  intact  (top),  25%  (mid)  and  50% 
(bottom)  tumor  volume  models  under  a  1  MPa  load  applied 
at  a  loading  rate  of  100  MPa/s. 
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Figure  7:  Three  dimensional  finite  element  mesh  of  a 
metastatically  involved  vertebral  body  and  adjacent 
Intervertebral  discs. 
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